Effects of C, Cu and Be Substitutions in Superconducting MgB 2 
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Density functional calculations are used to investigate the effects of partial substitutional alloying 
of the B site in MgB2 with C and Be alone and combined with alloying of the Mg site with Cu. The 
effect of such substitutions on the electronic structure, electron phonon coupling and superconduc- 
tivity are discussed. We find that Be substitution for B is unfavorable for superconductivity as it 
leads to a softer lattice and weaker electron-phonon couplings. Replacement of Mg by Cu leads to 
an increase in the stiffness and doping level at the same time, while the carrier concentration can 
be controlled by partial replacement of B by C. We estimate that with full replacement of Mg by 
Cu and fractional substitution of B by C, T c values of 50K may be attainable. 
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The discovery of superconductivity with critical tem- 
perature T c — 39 K in hexagonal MgB2 has led to con- 
siderable interest in this material, both for practical ap=. 
plications and from a more fundamental point of view.EI 
In particular, although it is not the highest T c conven- 
tional (in the sense of s-wave with likely substantial 
electron-phonon coupling) superconductor ,0 its substan- 
tial T c combined with its chemistry, and other properties 
may make it particularly usefulJjti A variety of experi- 
mental measurements, including B isotopj^cffcct,Ej nc 
tron scatt«™g-phonon densities of^atesjffi tunnelling, 
transport ,01130 and specific heatjl3t3 have been done. 
Taken together, they yield a consistent clear picture of 
MgB2 as a conventional electron-phonon s-wave super- 
conductor. This is also consistent with most of the re- 
ported theoretical work. 

The underlying electronic structure and bonding were 
discussed by Ripplingeio before the discovery of su- 
perconductivity. This was based on density functional 
calculations using the linearized augmented planewave 
(LAPW) method and provided results consistent with 
an early calculation by Armstrong and Perkins x3 Subse- 
quently, several authors have extended this work present- 
ing the electronic structureJ.n_detail and discussing the 
origin of superconductivity.E3~E3 The electronic structure 
is dominated by the chemical bonding of the hexagonal B 
sheets in MgB2. Although the nominal electron count is 
the same as in graphite (Mg has its nominal 2+ charge), 
unlike graphite, the top of the B derived bonding band 
structure contains holes related to fractional occupation 
of bands derived from B p states directed perpendicu=. 
lar to the layers. As pointed out by An and Pickett ,E£l 
the negative charge of the B layers and the corresponding 
positive charge of the Mg layers leads to a lowering of the 
p z bands relative to their position in graphite. Detailed 
calculations of the phonon spectra and electron phonon 
coupling have been reported. These show strong coupling 
between the hole "doped" bonding bands and high fre- 
quency optical phonons associated with motions of the B 
atoms affecting the covalent bonds. Kong et al. reported 
detailed calculations of electron-phonon couplings over 
the whole Brillouin zone obtaining couplings consistent 



with the measured T c as well as the boron isotope effect, 
the specific heat enhancement, the reported gap values 
and transport data.tZI 

The picture that emerges is one where the high T c is 
due to strong electron phonon coupling associated with 
the hole doped metallic a bonding bands in the B sheets. 
This is due to the strong bonding of these sheets, which 
along with the light B mass is responsible for the high av- 
erage frequency of the strongly coupled phonons, setting 
the temperature scale. However, there are also signif- 
icant interactions with lower frequency modes and, 
may be needed to explain the upper critical field! 
Though some interesting questions remain, e.g., regard- 
ing possible a nharm onic effects for the lower frequency 
phonon modes,E9uJ the basic physics appears to be set- 
tled as regards the origin of superconductivity in MgB2. 
The crucial aspects seem to be: (1) band structure, par- 
ticularly the presence of hole doped a bands at the Fermi 
level, Ep, (2) strong electron phonon coupling associated 
with the strong covalent bonding nature of these bands 
and (3) high phonon frequencies associated again with 
the strong covalent bonds and the light B mass. 

The present work is to explore some possible substi- 
tutions and provide suggestions for experimental work 
aimed at finding related high temperature superconduct- 
ing phases. Along these lines, Medvedeva et al. have 
investigated a number of possible substitutions of Mg by 
monovalent, divalent and trivalent ions.E3 They focussed 
on the band structure, and particularly the density of 
states (DOS) and presence or absense of the a band at 
Ep They concluded that trivalent substitutions like Y, 
Al etc. are not favorable as they fill the hole doped (in 
MgB2) a bands, while certain monovalent substitutions 
for Mg may be favorable. This is consistent with the 
experimental obsepjation that Al substitution destroys 
superconductivity.EHl They also mentioned the possibility 
of vacancies or substitutions on the B sheets, but con- 
cluded that these are all unfavorable. 

Here we briefly reexamine the effect of Be substitu- 
tion on the B planes, including band structure, lattice 
stiffness and electron-phonon calculations. Although Be 
substitution lowers the electron count, wc find that it is 



1 



so detrimental to the bonding that superconductivity is 
suppressed both due to a strong decrease in the lattice 
stiffness and a drop in the electron-phonon coupling. We 
then report calculations investigating the effect of a com- 
bined partial substitution of C for B and Cu for Mg. The 
rationale for this combination is that C-B bonds are ex- 
pected to be very strong in this structure, and so this 
substitution may lead to a stiffening of the sheets rela- 
tive to MgB2, while the replacement of Mg by Cu may be 
expected to first of all compensate for the extra charge 
provided by the substitution in the plane, and secondly, 
perhaps maintain the hole doping of the a bonding band, 
present in MgB2 but absent in graphite. The basis of this 
latter conjecture is that the crystal field in materials can 
arise both from strict ionic (Coulomb) effects, as seem to 
dominate in MgB2, and hybridization, which is impor- 
tant e.g. in transition metal oxides. The hope then was 
that such effects with monovalent Cu might be sufficient 
to keep the a band at the Fermi energy. As discussed 
below we find that both these effects do occur and that 
for modest C incorportation and Cu replacement of Mg 
superconductivity should be enhanced. Based on rigid 
muffin tin approximation (RMTA) calculations, we esti- 
mate that T c may be near 50K. 

Our results were—obtained in the local density ap- 
proximation (LDA)l3 using the several potential LAPW 
method, including local orbitals,cSc3 and well converged 
basis sets and zone samplings. Band structures, lattice 
stiffnesses and electron phonon couplings were all deter- 
mined at the calculated lattice parameters as determined 
by energy minimization. Commonly in LDA calculations, 
lattice constants are slightly underestimated. For MgB2 
we get a=5.736 ao and c=6.522.ao, as compared to the 
experimental values a =5.826 ao and c=6.653 ao. This is 
a 5% volume compression. As shown in Fig. [I] this has a 
very small effect on the band structure near the Fermi en- 
ergy, Ep with a small effect on N(Ep). To proceed to the 
alloys and discuss superconductivity, we make three more 
significant approximations. First of all, we employ the 
virtual crystal approximation (VCA) to account for par- 
tial substitutions on the B sheets. A partial justification 
for this is provided by the strong covalency and corre- 
sponding large bandwidths, which may limit the amount 
of scattering due to potential disorder in the alloy (this 
is the same effect that allows alloys, like Al^Gai-^As, 
to have high enough mobilities to be useful in semicon- 
ductor technology). We tested the VCA by comparing 
with ordered cells at the compositions MgBeB and CuBC 
and found some quantitative differences (discussed be- 
low), but the key features of the band shapes, velocities 
and the position of the a bonding band relative to Ep 
were little changed. The second approximation we made 
was to characterize the lattice stiffness that combined 
with the mass gives average phonon frequencies by the 
calculated tensile stiffness of the B sheets, i.e. d 2 E/da 2 . 
This number and the average mass of the atoms in the 
B sheets was used to scale the average frequencies as cal- 
culated by Kong et alxJi Considering that the dominant 



phonons determining T c are the B modes, we think that 
this is a reasonable approximation. Finally, we use— the 
RMTA to characterize the electron phonon couplingESEZl 

We used well converged DOS calculations based on first 
principles eigenvalues at over 2000 k points in the irre- 
ducible zone to get N(Ep) and the corresponding angu- 
lar momentum components. We used the self-consistent 
LAPW potentials at different concentrations to calculate 
the corresponding phase shifts and free scatterer DOS at 
Ep. Sphere radii of 1.5 ao were used for B, Be and C for 
the RMTA calculations. The above quantities were used 
in the Gaspari-Gyorffy formula to compute the Hopfield 
parameters, ij, for each site. Negligible coupling is found 
on the Mg site, as expected, but not Cu, e.g. 0.78 eV/A 2 
on Cu for CuB2. However, in Table | we give the values 
of rj for B only, since that is the dominant contribution 
controlling superconductivity and we do not include any 
Cu contribution in the calculation of A or T c (Cu will 
have little involvement in the high frequency phonons 
associated with the B sheets). For the electron-phonon 
coupling we used the usual expression A = rj/< Mlo 2 >. 
For the denominator we-psed the average frequency cal- 
culated by Kong et alS-i for MgB2 and scaled it using 
the tensile stiffness of the B sheets for various concentra- 
tions, i.e. d 2 E/da 2 . The RMTA is not generally as well 
justified in sp metals as in transition metals and can con- 
siderably underestimate the deformation potentials when 
strong sp covalent bonding is present as it is here. Fur- 
ther, the RMTA neglects some differences between dif- 
ferent bands, which may be significant here. In any case, 
Kortus et al. did use it for MgB2 to characterize electron 
phonon couplings .113 Comparing with the direct calcula- 
tions of Kong et al. for MgB2 we find, not unexpectedly, 
that the values of 77 we obtain with our non-overlapping 
B spheres are too small, roughly by a factor of three. 
Here we use the RMTA to elucidate trends, while ac- 
knowledging its limitations. So we adopt the heuristic of 
scaling the calculated values of 77 by three in calculating A 
and T c . Using the average phonon frequencies quoted by 
Kong et al. this heuristic closely reproduces their values 
of A and T c for MgB2. We used the McMillan equation to 
roughly estimate T c , setting the Coulomb pseudopoten- 
tial n* =0.1. We emphasize here that we are not aiming 
at an accurate determination of the value of T c but ex- 
ploring the trends upon substituting Mg by Cu and B 
by Be or C. Because of our use of the RMTA our values 
should be taken as indicating trends not as quantitative 
predictions of T c . 

The band structure for the starting compound, MgB2, 
(Fig. [l])|-jS|— practically identical to those obtained 
previously|l3il2l'E3 showing a bonding states at Ep. Re- 
sults for the structural and electronic properties relevant 
to superconductivity are given in Table Q, while band 
structures are shown in Fig. || and Fig. || for MgBe a; B2_ a ; 
and CuB2-a;C 2 ;, respectively. In the table the values of 
77 are the bare values as given by the RMTA, while the 
heuristically scaled values have been used for the values 
of A and T c . An ordered band structure for CuBC is 
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shown in Fig. [|. This shows some differences from the 
corresponding VCA calculation, most notably, near Ep, 
a splitting at the H point involving p z bands. However, 
the general structure of bands near Ep and the position 
of the a band is quite similar. 

Substitution of Be into the B sheets lowers the electron 
count, though as seen in the band structure, not in a rigid 
band way. The hole concentration in the a bonding band 
does increase, but the bonds are strongly weakened. This 
is reflected in the band widths and lattice stiffness. The 
result (Table is a rapid increase in the a lattice parame- 
ter, softening of the lattice and a decrease in the electron 
phonon coupling. Thus Be substitution is detrimental to 
superconductivity and is not further discussed. 

The Cu substituted material is more interesting. In 
an ionic model, replacement of Mg by monovalent Cu 
should lower the electron count in the B sheets by one 
per formula unit. So perhaps a situation similar to that 
found with 50% Be substitution may be expected. This 
is not the case. The Cu is monovalent as in the ionic 
model. The five narrow Cu d bands are in the valence 
region between -4 and -3 eV relative to Ep. Comparing 
the top panel of Figs, [l] and || and the structural infor- 
mation in Table [| one sees that the bonding in the sheets 
is strengthened by the addition of Cu, even though the 
hole count in the a band is increased (note the relative 
positions of the band maxima on the T — A line). Hy- 
bridization with Cu leads to quantitative changes in the 
band structure, affecting mostly the p z states, but there 
are also weak effects on the a bands near Ep, e.g. the 
reversed dispersion on the T — A line presumably due to 
d—pn interactions. The net effects of Cu substitution - 
stiffened lattice, increased N(Ep) and higher hole con- 
centration - are favorable for superconductivity. Related 
to this, there is a recent unconfirmed report of an ep-. 
hancement of T c with partial replacement of Mg by Cu.Ea 
Partial substitution of C for B in CuB2 has two effects 
- a stiffening of the lattice reflecting the strong bonding 
of C and B (favorable for T c and a reduction in the hole 
doping of the a band and in N(Ep) (unfavorable for T c ). 
For low C concentrations, the first effect dominates, lead- 
ing to an increase in the estimated T c , but beyond 25% 
C substitution (here the van Hove from the a band at A 
crosses Ep) the second dominates and T c falls. 

In summary, band structure and RMTA calculations 
indicate that the T c of MgB2 can be increased, perhaps 
to 50K by substitution of Cu for Mg and low C substitu- 
tions, around 25%, in the B sheets. We emphasize that to 
our knowledge these are not known equilibrium phases. 
However, the strong binding of the sheets suggests syn- 
thesis may be possible, e.g. starting with MgB2 (n.b. Mg 
leaves the compound at high temperature). 

This work was supported by ONR and the the ASC 
supercomputer center. We are grateful for discussions 
with I.I. Mazin, W.E. Pickett and K. Schwarz. 
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28 Note that N(E F ) in Table § is at the LDA lattice parame- 
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TABLE I. Calculated properties of MgBe^E^-a; and 
CuB2- a; C I as obtained in the VCA. S is the relative stiff- 
ness of the B sheets characterized by d 2 E/da 2 (S(MgB2) = l), 
and the other symbols have their usual meanings. Lattice pa- 
rameters are in ao, N(Ef) in eV _1 and T c in K. The unsealed 
7] (in eV/A 2 ) is for the alloyed B site only. Blank entries for 
a, c and S indicate that these were obtained by interpolation 
rather than direct computation for those concentrations. 
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FIG. 1. LDA band structure of MgB 2 with the ex- 
perimental structure (solid) and calculated LDA lattice 
parameters (dashed). The zero is at Ef- 





FIG. 3. LDA virtual crystal band structures of 
Cvfi$2-xCx for x = (top), x = 0.5 (middle) and x = 1 
(bottom). The lattice parameters are the calculated re- 
laxed values. Note the vertical scale of the lower panel. 
The horizontal reference at denotes Ef- 



FIG. 2. LDA virtual crystal band structure of 
MgBeajF^-z for x = 1. The lattice parameters is the 
calculated relaxed value. The horizontal reference at 
denotes Ef- 





FIG. 4. LDA ordered band structure for CuBC note 
the splittings relative to the virtual crystal band structure 
in the bottom panel of Fig. |j| 
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